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Abstract 

A quantization p root'd u re for the electromagnetic field in a rectangular cavity with perfect 
conductor walls is presented, where a decomposition formula of the field plays an essential 
role. All vector mode functions are obtained by using the decomposition. After expand- 
ing the field in terms of the vector mode functions, we get the quantized electromagnetic 
Hamiltonian. 


1 Introduction 

Recently we have earritxi out the field quantization in several rectangular cavities using the vector 
mode functions [1]. The vector mode functions have been obtained with the help of an orthogonal 
matrix. However, the procedure developed there has not been applicable to other cavil, it vs in a 
straightforward manner. 

To overcome tint above difficulty, we have presented another quantization scheme for the field 
in a circular cylindrical cavity [2]. All vector mode functions have been obtained by using a 
decomposition formula derive! from Maxwell's equations directly. This method is more general 
than before, because if is also applicable to rect angular and spherical cavities. 

In this paper, we applied the above method to a rectangular cavity with perfect conducting 
walls. Then spontaneous emission of an atom inside the cavity is calculated. 


2 Decomposition of Electromagnetic Fields 

In this section, we derive the decomposition formula for the field in the Cartesian coordinates. We 
shall use this nsailt in performing the field quantization in a rectangular cavity in Sec. 3. 

Maxwell’s (equations for the electric field E and the magnetic field B in free space are given by 
V ■ E 0 , V • B 0 , and 


V x E f dfcB 0, V x B - J^E -- 0, (1) 

where c is the velocity of light in free space and d/dt. 

The electromagnetic field E and B can be written in the Cartesian coordinates (x, y, z) as 
E E r I E 2 , B B r I B 2 , where E r e x E x \ e y E y is the transverse component of the field 


223 


and E z e z E z . Mere e x , e y , and e 2 are the unit vectors in the x, y, and 2 directions, respectively. 
For simplicity, the derivative with respect to, for example, x is described as d/dx -- d x . 

Since V 2 x E 2 0, the first equation in (1) gives 

dt B 2 — Vt x E T , (\Br Vt x E 2 V 2 x E^. (2) 


Similarly, the second equation in (1) leads to 

—dtE z — Vr x B t , -^Et Vt x B 2 ( V 2 x Bp. (3) 

c 2 c 2 

Equations (2) and (3) give 

A t E - -V x V x E 2 + V x dtB z , A t B = — x #E 2 - V x V x B z . (4) 

& 

To rewrite Eq. (4), we must, decompose the components E z and B z into two parts. Suppose 
that the field is in a finite region. Let us expand E z and B z in terms of a certain complete system 
of functions with mode s: 


E z (r,t) ^(/? 2a (r,0 +c.c.), B z (r,t) = ^(^(r,*) + c.c.), (5) 

a a 

where E za {r,t) £ 2S (r)e~ ,w * lt and B za (r,t) - B 2a (r)e _, "* 2< . Here uw (u M > 0, a - 1,2) is 
determined by using given boundary conditions. Since E z and B z satisfy the wave equation, the 
components E M and B za satisfy the Helmholts equations: 

A E„ -kl.Ev,, A B^ - -l&B za ( r), (6) 

where k air u^/c 2 . We assume that the components satisfy 

%E za - -h 2 al E za , <? Z B „ = -h 2 a2 B za , (7) 


where is determined by the boundary conditions. Then we have two dimensional Helmholtz 
equations: 

A T £ 2a (r) -g 2 al E za ( r), A T B za (r) - -&B„( r), (8) 

where g 2 a<r k 2 , - h 2 M . 

Here we define two functions F„ from E za and B za with ^ 0 as 

K(r,l) - E |ft.(r,l) +c.c.j = E IF-Me-”-* +c.c.], (9) 

9io^° 9ia^° 

where 

F al = E„/£ u Fa = B m l&* (<& / 0 ). ( 10 ) 

The functions F„ and their components F M satisfy the same equations as E z and E^, respectively. 
The component F a is a solution of the Poisson equation. On the other hand, if there is a component 
E za or B za with 0, Eq. (8) reduces to two dimensional Laplace equation. 
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The functions A r E n satisfy 


— A t Fi E .0*i 

e F 2a , 

—At ^2 5^ 9a 2^*2 ^zs- 

(ii) 

9»l A) 




I Iere define E )z and B( )z as 




Eo z E z \ At Ei 

E 

Boz — B z + At/^2 = B zs , 

(12) 


sii=° 

£ 2^° 


which satisfy AtEq z 0 and AtBq z 

0. Then 

we have useful formulas for E z and B z \ 



E z - A/-/') ( Eq z , B z — —A f Bq z . (13) 

Using Eq. (13) and defining F„asF„ e 2 F„, we can rewrite Eq. (4) as 

A t (e - V x V x Fi -t Vx &F 2 ) = -V x V x Eoz + V x dtB 0z , (14) 

A^B-^Vx^F.-VxVx F 2 ^ -^V x ^E 02 - V x V x B 0z , (15) 

where E<> 2 e z E 0z and B 0z e z B 0z . Define E 0 and B 0 as the quantities in the parentheses at 
the left hand side in Eqs. (14) and (15), respectively. The results of this section is summarized in 
the following theorem. 

Theorem 1: If the components E za and B za satisfy Eq. (7), the field can be decomposed into 
three components as follows: 

E VxVxFi-Vx dfF 2 4 Eo, B -^V x $F, 4 V x V x F 2 + B 0 , (16) 

where E 0 and B 0 satisfy 

A t E 0 -V x V x E 0 z 1 V x dtB 0z , A r B 0 = -^V x dtEo z - V x V x Bq 2 (17) 

Theorem 1 plays a central role in performing the field quantization in this paper. It is worth 
emphasizing that, each term in Eq. (16) is a solution to Maxwell’s equations. 

3 Vector Mode Functions and Field Quantization 

The cavity we treat here is enclosed by rectangular walls having sides L x , L yy and L z in the x, y, 
and z directions, respectively: 0 < x < L x , 0 < y < L y , and 0 < z < L z . We assume that the 
cavity has perfectly conducting walls. The tangential component of the electric field E| tan and the 
normal component of the magnetic field B| norm must accordingly vanish at the boundaries of the 
cavity. 
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Tin* above boundary condition reduces to that for the 2 components 

E z 0, () X B Z 0, (z - 0, L x ), (18) 

E z 0, d y B z 0, (y 0, L y ), (19) 

B z 0, d z E z 0, ( z = 0, L z ). (20) 

The solution to the Helmholts equation (6) for the components Ez, and B za under the above 
boundary conditions is given by 

E za (r,t) C 9 i(t) sm(£nx/L x ) sin(mny/ L y ) cos(mrz/ L z ), 

B za (r,t) C a2 (t) cos(£nx/ L x ) coe(mnz/ L z ) s\n(nnz/ L z ), (21) 


where the mode index is s (£, m, n) (£, m, n — 0, dtl, ±2, • ■ •). 
fVorn the solution (2 1 ) we have 

air = g 2 a {tn/L x y 4 (rwr/L,) 2 , 

kl = k 2 (tir/L x ) 2 + (mTr/4) 2 + (nir/ L z ) 2 . (22) 


Consequently, we can use Theorem 1 in the preceding section. Although it follows from Eq. (22) 
that > 0, we can prove that g\ a > 0, which results from the following lemma. We omit its 
proof. 

Lemma: Eq, Bq 0 and gl a >0. As a result, the term with i m = 0 in Eq. (21) cannot be 
used. 


Let us next obtain the functions F„, whose definitions are given in Eqs. (9) and (10). That is, 
the functions are given by 


F„,(r,i) 


F a2 (r, t) 


E za (r,t) . j hu) a 

2e 0 


B. 


/ 


hu) a 

2e 0 


a ai (t)ip a i(r). 


a a 2 (t)if} a 2 (r), 


(23) 


where u) aa = u>, and we have introduced 0 ^( 1 ) — a a<r ( 0)e *“'** and ip*, given by 


ip a i (r, t ) — c a i sin(tnx/ L x ) sm(mny / L y ) cos(nirz/ L z ), (24) 

xp a2 (r,t) c a2 cos(£nx/ L x ) cos(m7ry/ L y ) sm(nirz/ L z ), (25) 


with c a i [8/(V k 2 g 2 )\ l/2 and c a2 - [8 / {V ulg 2 )\ l/2 . The functions xp a „ has the orthonormality 
property 

(26) 

where f c dr / cavjty dxdydz and V is the cavity volume. Here the quantity cosrr in Eqs. (24) and 
(25) must be changed to 1/ \/2 when x — 0. 

Substituting the functions F„ in Eq. (23) into E in Eq. (16), we find 


E(r ,0 *St/^[«^(0*w( r )-«lr(0 , C( r )]» 

M V 


(27) 
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where the vector mode functions u*, are given by 


u«i V x V x e z ip„i, u a2 iuJ 9 V x e z i/> a2 . (28) 

'I’he vector mode functions satisfy V • 0, (A t k a<T ) u,„ 0, and at the boundaries u a<7 | tan - 

0, V x u m | norm 0. They also satisfy the orthonormality property needed for quantization: 

J dr u^(r) ■ u a ' a > (r) - 6 M > (29) 

To get the quantized field, the functions a a „{t) are regarded as annihilation operators satisfying 
the commutation relation [a*r(t), “ bamboo’- Then we get the following theorem. 

Theorem 2: The quantized field and the Hamiltonian are given by 

E(r,() u„(r)-«l.(l)u^(r)], (30) 

B(r, 0 \l^~ [<w(*)V x u^r) T «L(t)V x u^(r)] , (31) 

//fi ^ ~ fad a 1 ~ ^ 2) ' (^^) 


4 Spontaneous Emission 

As an application, we consider the transition rates of an atom in the cavity, using the dipole 
approximation. The Hamiltonian is given by H — Ha t Hu 1 Hi, where Ha is the free Hamiltonian 
for the atom, Hr for the field which is given in Eq. (32), and Hi = eD ■ E(R) (— eD: the total 
electric dipole moment of the atom; R ( X , Y, Z): the position of the atom). 

At t --- 0, the atom is in an energy state |i 0 > (with energy E^) and the field is in the vacuum. 
Then the probability per second of finding the atom in a state (with energy E t ) at sufficiently 
large time t is given by 


w 


e 2 7r 

£oh 


£|iUR)-<i 0 |D|»>| 


2 cJ a sin(cJ, — u}p)t 
ir(uj a — u>o) 


(33) 


where hu ) 0 Ei 0 — E t . 

Let us take the average of the coordinates Y and Z and take L z —* 00 . The transition rates 
w y and w z vanish ( Wi indicates the rate where the dipole moment is along with the i direction). 
The rate w x is given by 

wjm 0 AD 1 " - ">73). (34) 

^sy m 

where & u 0 Li/cK, £ x -- 1/2, Wq - e 2 | < iojOji > \*u)QH$mh£Q(?), 0(x) = 1 (x > 0), and 
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0(x) 0 (x < 0), which is shown in FIG. 1. 



FIG. 1. Transition rate w x for L z — > oo, where the dipole moment is along with 
the x direction. 

Setting here u 0 10 13 Hz, we have L x = 4.7 x 10~ 2 mm, so that the cavity is quite narrow. 
Also, FIG. 1 shows that the transition is forbidden when < 1, i.e., L y < 9.4 x 10 -2 mm., where 
the cavity is a thin tube in this case. 

5 Conclusions 

The quantization for the field in the cavity has been performed as follows: obtain the decom- 
position formula (16) in the Cartesian coordinates; solve the Helmholts equations (6) for the 
components E za and B za under the boundary conditions; determine the functions F v , substitute 
them into the decomposition formula (16), and obtain the vector mode functions satisfying the 
orthonormality property (29); then we arrive at the quantized field and Hamiltonian. In the whole 
process of quantization, the decomposition formula in Theorem 1 plays an important role. 
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